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Abstract 

We find the analytic expression of trp"(I/) for a massive boson field in 1 + 1 dimensions, where p{L) 
is the reduced density matrix corresponding to an interval of length L. This is given exactly (except 
for a non universal factor) in terms of a finite sum of solutions of non linear differential equations of 
the Painleve V type. Our method is a generalization of one introduced by Myers and is based on the 
explicit calculation of quantities related to the Green function on a plane, where boundary conditions 
are imposed on a finite cut. It is shown that the associated partition function is related to correlators 
of exponential operators in the Sine-Gordon model in agreement with a result by Delfino et al. We 
also compute the short and long distance leading terms of the entanglement entropy. We find that 
the bosonic entropic c-function interpolates between the Dirac and Majorana fermion ones given in 
a previous paper. Finally, we study some universal terms for the entanglement entropy in arbitrary 
dimensions which, in the case of free fields, can be expressed in terms of the two dimensional entropy 
functions. 



1 Introduction 

The trace of the vacuum state projector over the degrees of freedom lying outside a given set A gives 
place to a local density matrix p^j, which is generally mixed. The associated entropy >5'(^) is called 
geometric or entanglement entropy. In recent years there has been much interest on the properties of the 
local reduced density matrices. This is partially motivated by ideas related to black hole physics (i.e. see 
^ 121 El E]) I developments in quantum information theory, and the density matrix renormalization group 
methods in two dimensions (within the area of condensed matter physics see for example 0E1)- Besides, 
the entanglement entropy as well as other measures of information for the reduced density matrices of 
the vacuum state are in its own right interesting and relevant quantities in the context of quantum field 
theory. For example, they show a different structure of divergences and a nice interplay between geometry 
and ultraviolet behavior. Among the results in this sense we can mention the existence of an entropic 
c-theorem [7] and the recent discovery of a topological entanglement entropy |S] . 

In a previous paper we studied universal quantities derived from the local density matrix pA, for a 
Dirac field in two dimensions j^. In the present work we find the corresponding results for a massive 
scalar field. Specifically, we shall consider the a entropies 

^o(^) = 7^1ogtr(pl), (1) 

and the entanglement entropy 

S{A) = -tripA log pa) = lim 5„ (A) , (2) 

where A is a single interval of length L. As these quantities are not well defined in the continuum limit, 
we work instead with its associated dimensionless and universal functions Ca{L) and c(L) 

c„(L) ^ L^^^ , c{L) ^ ci(i) . (3) 
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The entropic c-function c(L) is always positive and decreasing, and plays the role of Zamolodchikov's 
c- function in the entanglement entropy c-theorem 

The traces trp^ involved in with a = n S Z, can be represented by a functional integral on an 
n-sheeted surface with conical singularities located at the boundary points of the set A [3 Ej. Here, 
calling u and v to the end points of the segment, the replicated space is obtained considering n copies of 
the plane cut along the interval {u,v), sewing together the upper side of the cut {u,v)^^^^ with the lower 
one for the different copies fc = 1, and where the copy n + 1 coincides with the first one. 

The trace of p" is then given by the functional integral Z[n] for the field in this manifold, 

trp" = 4^ . (4) 

We consider a complex scalar field and following [3, we map this problem into the one of a single 
cut plane with n decoupled fields <I>fc(a:), k = 0, — 1. These are multivalued, since when encircling 
the points m or z; in anti-clockwise sense they get multiplied by e'™^'^ or e~'"^'^ respectively, with k = 
0, n — 1. Then 

n-l 

logZ[n]^Y.^ogZk/n- (5) 

k=0 

where Za, with a G [0, 1], is the partition function of a complex scalar field in a plane with boundary 
condition 

= e'2.''^out(x) (6) 

along the cut {u,v). 

The free energy and the Green function are related by the identity 

d„,2\ogZa = - J dr^G{r,f). (7) 

Our strategy is to find the right hand side of this equation. The information on the Green function 
on a cut plane which is relevant to the present problem is obtained following a generalized version of a 
method introduced in to deal with Neumann and Dirichlet boundary conditions on the finite cut. It 
essentially consists in exploiting the rotation and translation symmetry of the Hclmholtz equation, even 
in the presence of the cut boundary conditions, by analyzing the behaviour at the singular points. In 
Section 2 we adapt this method to our particular boundary conditions and find an exact expression for 
dhdm'^ log Za in terms of the solution of a second order non linear differential equation of the Painleve 
V type. This equation allow us to identify the partition function with the inverse of a correlator of 
exponential operators in the Sine-Gordon (SG) model (a related mapping for correlators of boson disorder 
operators was given in jllp. 

This is reminiscent to the fermionic case, where there is also an expression for Ca in terms of SG 
correlators. This is shown in a previous paper by a very different method involving dualization jHj. 
The relation between the fermionic and bosonic cases is analyzed in detail in Section 3, where we first 
present the results for the bosonic a entropies. In this Section we also study the short and long distance 
expansions for the entanglement entropy, and we show that the entropic c-function for a real scalar 
interpolates between the ones corresponding to a Dirac and Majorana fields. In the Section 4 we present 
some results concerning universal terms for the geometric entropy in arbitrary dimensions. These can 
be exactly computed for free fields in terms of the two dimensional entropic c- functions. The Section 5 
contains the conclusions. 

2 The Green function 

In this section we study the Green function G{f,r') using a method introduced by Myers in jlOj. The 
following three requirements uniquely define G{f,r'): 

a.- It satisfies the Helmholtz equation 

{-Ap + m^)G{r,f')=S{r-^). (8) 
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b.- The boundary conditions are (they also hold for the Green function derivatives) 

lim G((x,e),r ) = e^^'^Mim G'((a;, -e)^ ) ioi x e [L2, Li] , (9) 

e— »0+ e— »0+ 

]im G{f,/) = 0. (10) 

\r\ — 'oc 

Here the interval [L2, Li] on the x axis is the cut location. From now on we choose a € [0, 1]. 
c- G{r,r') is bounded everywhere (including the cut) except at f = r*. 

We will write the Green function as G{z, z') as a shortcut of G'(z, z, z', z', Li, L2), where z and z are the 
complex coordinates x + iy and x ~ iy. It is Hermitian 

G{z,z')^G{z\z)* , (11) 

and the symmetry of the problem also gives 

G{z,z'y ^G{z,z'). (12) 

The reflexion operation 

R{x,y)^{Li+L2-x,y) (13) 
leaves the Helmholtz equation, the cut, and boundary conditions invariant. Thus we have 

G{z,z')=G{Rz,Rz'). (14) 

Due to the boundary conditions, near the end points of the interval [L2, Li] the Green function must 
have branch cut singularities. The requirement that the function must remain bounded on the cut and 
the equation JHJ) imply that the most singular terms of G{z, z') for z near Li (and fixed z') have to be of 
the form 

G(z, z') ~ (z - LirSiiz') + (z - ii)i-'^52(z') . (15) 

We have written explicitly only the terms with powers of z — Li with exponent smaller than one. These 
are the ones of interest in the sequel. Note that the contributions at this order must be analytic or 
anti- analytic in z in order to cancel the Laplacian term in (|SJ|. 

The most singular contributions to G(z, z') for z Li follow from the derivative of H15|l 

5l,G(z, z') ^ -a(z - ii)'^"i5i(z') + (a - l)(z - Li)-^S2iz') . (16) 

The function dLiG{z, z') satisfies the equation ijSJ and the boundary conditions, and it is not singular 
at z — > z'. It has only one singularity located at ii whose expression is given by (|16|l . Therefore the 
combination 

9l,G(z,z')(5i(z")52(z"') - 5i(z"')52(z")) +aL,G(z,z")(5i(z"')52(z') - 5i(z')52(z"')) 
+9^, G(z, z'") (5i {z')S2 (z") - 5i {z")S2 (z')) (17) 

is nonsingular everywhere, and the uniqueness of the solution implies that it must be identically zero. 
Thus, 

aL,G(z,z') = ^^i(z)5i(z') +F2(z)52(z') = F^{z)S{z') (18) 

for certain functions -F'i(z) and F2{z). Note that in the last equation we have written the sum over two 
terms in vector notation. From the Hermiticity relation we have 

S'(z) = A-iF(z)*, (19) 

where A is a Hermitian matrix. Finally we get the following fundamental relation 

9l,G(z, z') = F'^{z)S{z') = Siz^ASiz') . (20) 



3 



The function S{z) satisfies the homogeneous HelmhoUz equation with the same boundary conditions 
as G{z',z) with z' fixed, but it is unbounded around Li. fn fact, it follows from Hlf)|l and 12()|l that it 
has singular terms proportional to (z — and (2 — However, these disappear if we derive 

with respect to ^2- On the other hand S{z) behaves around L2 as an ordinary wave, that is, it vanishes 
proportionally to {L2 ~ z)" and {Li — z^^"" . Then dh^iS^z) has singular terms around L2 which are 
proportional to (L2 — z)"^^ and (L2 — z)^". These have the same form as the singular terms of S{Rz) 
for z near L2 once one takes into account that 

z-Li=L2-Rz. (21) 

In consequence, there exist a matrix 7 such that d^^ S{z) — 7 S{Rz) is bounded, and thus identically zero, 

dL,S{z)^^S{Rz), (22) 
where 7 is a function of L = Li — L2. Translational invariance then implies 

dL,SiRz,Li,L2) = idi+d2){S)iRz,Li,L2) = -d3{S){Rz, L^, L2) = -7^(^) • (23) 

Using (pn|l . and (|^ to calculate dLidL2G{z, z') — dL2dLiG{z, z') = we get that the matrix A must 
be a constant, dAjd^i — dAjd^^ = 0, and it also holds 

7!' = A-iA-^ . (24) 

The constant matrix A can then be obtained using the solution for the half infinity cut, which can be 
computed by standard methods. In polar coordinates the Green function with L2 = — cxd, Li = is 

G(r, 0, r\e')^— J2 e'("+'^)''e-'("+'^)''7,„^„, (mr<)if|„+,| (mr>) , (25) 

n— — oo 

where / and K are the standard modified Bessel functions and r> and r< are the maximum and minimum 
between r and r' . From here we get 

A = -47r(l - a)acri , (26) 
where ai is the Pauli matrix. The equation H16|) leads to the behavior 

for z in the vicinity of Li. 

2.1 Equations for S{z) from translation and rotation symmetries 

In order to use eq. H2UI) to compute the partition function we need more information on S{z). With this 
aim we exploit the symmetries of the problem as follows. Translations in the y direction commute with 
the Helmholtz equation and thus the function {dy + dy')G{z, z') is non singular at z = z' and satisfies 
the equation © plus the boundary conditions. However, due to the dyG{z, z') term, it is singular at Li 
and L2 as a function of z. Near Li we have from (|15|) . H19|l and H27|l 

dyG{z, z') - ^a(z - Li)^-^Siiz') - i{l - a)(z - Li)-''S2iz') - -iFi{z)Siiz') + iF2{z)S2{z') , (28) 

and near L2 we have by the symmetry (I14II 

dyG{z,z') - ~iFi{Rz)Si{Rz') + iF2{Rz)S2{Rz'). (29) 

Then, by subtracting these terms from [dy + dyi)G{z, z') we obtain a singularity free function which 
respects the boundary conditions and satisfies the Helmholtz equation. In consequence it is the null 
function, 

dyGiz, z') + dy>G{z, z') = -i(F^(z)CT35(z') + F^ {Rz)asS{Rz')) . (30) 
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Taking the derivative with respect to Li we arrive at 

dyF^iz)Siz') + F^{z)dy,Siz') = i{F^iRz){-f,a3}S{z')+ (31) 
F^{z) {7, ^3} S{Rz') + d^F^(z)asS{z') + F^ {z)a^d^,, S{z')) , 

with {7, tJa} = 7(T3 + (T37. This equation has to be satisfied in z and z' independently. In particular 
for z ^ Li the coefficients of the singular terms on both sides of H31() must be equal. From H19() and 
(I27|l the most singular terms in dyF^{z) — idxF^ {z)a^ (which would be proportional to (z — LiY~^ and 
[z — Li)~"~^) exactly cancel. But this combination of derivatives should also have a contribution to the 
next orders, (z — Li)°-~^ and (z — for z near Li. Then we can write near this point 

a^F^(z) - id^F'^{z)cj^ ^ F^{z)^ , (32) 

where ^ is a matrix which depends on L. Thus, isolating the singular terms of 13111 for z near Li, we 
have 

dyS{z) = I {7, ^3} S{Rz) + ia^dxS{z) - ^S{z) . (33) 
Replacing this back into (|31|l and using {A, 0-3} — 0, {7, (73}^ A + A {7, 0-3} — 0, we obtain 

+ AC = . (34) 

We can also write the reflected equation for (|33|l 

dyS{Rz) = i {7, 0-3} S{z) - iasd^SiRz) - (S{Rz) . (35) 
Taking the derivative of H33|l with respect to y and using H33|l . (|35|l and the Helmholtz equation we get 
(to2 -,_ ~ e) S{z) = -t {e, as} d^S{z) - I {{7, (73} , a S{Rz) . (36) 

The structure of singularities in this equation is different for the three terms, which have to cancel 
independently. This leads to 

{C^a} = 0, (37) 

{{7,'^3},^} = 0, (38) 

' + {7,^3}'-^') = 0. (39) 



The same line of reasoning followed above to treat the translation symmetry in the y direction can 
be used to exploit the rotational symmetry of the Helmholtz equation. Considering rotations around L2 
with operator {de — L2dy) where dg = xdy — ydx, we have that {{dg — L2dy) + {dg' — L2dy')) G{z, z') is 
not singular for z — z' . Analyzing the singular terms we get the equation 

{{dg - L2dy) + [dg, - L2dy.)) G{z, z') = ~iL F^{z)a3S{z') . (40) 

Taking the derivative with respect to L2, and following the same steps as before we arrive at 

{dg - L2dy) S{Rz) = -iqS{Rz) + iL-iu^Siz) , (41) 

where 

1-0^ m 

The reflected equation follows from replacing z — > Rz and {dg — L2dy) — {dg — Lidy), 

{dg - Lidy) S{z) = iqS{z) - iL-fa3S{Rz) . (43) 

We derive (|41|) with respect to Li, and use H33|l to bring it to the form of the equation H43|l . From 
the comparison we obtain the following matrix differential equation 

^^^(Lj-^^(73+j-\-r + q + a3j. (44) 
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2.2 Parametrization 

The algebraic equations H24|l . (|34|l . (|37I - I39|I for the matrices are solved using the parametrization 

where u, b, c, f3i, and (32 are real functions of t — mL and 

u^ + l^ (3iP2 ■ (46) 

The differential equation H44|) writes in components 

tbu' — tb'u — ub 

t[u'^ — be) 
tcu' — tc'u — uc 
t{u^ - be) 

= {1 -2a)u'^ +tuu' -tc'b + 2{a-l)bc, (49) 
= (1 - 2a)u'^ - tuu' + tcb' + 2abc. (50) 

Defining h ^ be the last two equations give 

h+^h' = tuu'. (51) 

Using 1I49() and H50|l we eliminate e'b and eb' in terms of u and h. Then H46() gives /i in terms of u and 
using H51|l it follows that 



/5i = ^t:^^^::^' (47) 
/32 = ^77^^::^' (48) 



u" + i«'--^^'2-,i(l + «2)_^i^^^0. (52) 



This nonlinear ordinary differential equation can be transformed to take the form of a Painleve V equation 
\12\ . Taking into account that u must decay exponentially fast at infinity we also have 

■ ,,2 t'u''-umi-2a)'+tHl + u')) 

dssuis) = . (53) 

2.3 Integrated quantities 

The equations Q and (|20|l give for the partition function 



dLd,n^ lO! 



where 



gZa = ~ J S'U 5 = Stto {l-a)H, (54) 

H{L) = j dxdy Sl{z)S2{z) . (55) 

To find this quantity, we use the information obtained in the preceding Section. We first define the 
following auxiliary integrals 

(56) 
(57) 
(58) 
(59) 
(60) 





j dxdy Sl{z)Si{Rz 


B2{L) = 


j dxdy S2{z)S2{Rz 


Bi2{L) = 


j dxdy S2{z)Si{Rz 


X,{L) = 


j dxdy 5*(z)5i(z), 


X2{L) = 


j dxdy S2{z)S2{z) . 
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These are convergent. They are also real since the relation (|12(l implies that S*{z) — Si{z). Combining 
(jH^II and H,S5|1 we obtain 

~i [dy - id.^) (SiS;) = mu{S^S; + SiS;^) - m {(iiS2S* + P^SiSl) . (61) 

Integrating this equation in the plane it is 

-i I dxdy {dy - id^) [SiSD = = 2muBi2 - m {P1X2 + /32^i) , (62) 

J ona (1 — a) 

where we have taken into account the pole of S1S2 at Li 

SiS; ^ ((4^)' a (1 - a)) (z - L^y' , (63) 

and the fact that 5*1 6*2 and its derivatives are continuous along the cut. Similarly, several other equations 
for the integrated variables can be constructed. From the angular equations (|41|l and (|43|) we get 

i J dxdy (de - L^dy) {SiSD = tuB^2 - \{hB2 + cB^) + H . (64) 

Using the eqs. and H35|) we obtain an expression for J dxdy S2*dySi in terms of Bi and B2, and a 
different expression for the same quantity arise using the angular equations l|41|) and H43|l . Combining 
them we get 

0=^(1- 2a)Bi2 ~ cXi + bX2 - P1B2 + P2B1 . (65) 

We still need three more equations involving the integral variables to completely determine them. A 
new one results as follows. It is easy to show that 

Jsixd^S, = ^^ = JsiydySi, (66) 

S;xd,S2 = -^ = Js;ydyS2. (67) 
Using the angular equations we have 







SiyS; = -i J S^xSi , (68) 

and from the eqs. and using and ^7^, 

J SlydxSi = -J SlxdySi ^ + Pirn j SlyS2, (69) 

j S;yd.,S2 = -J S;xdyS2 = -t^ + j S;xSi . (70) 

Now, inserting back these relations into the angular equations and taking into account Ht)8|l we finally 
have 

-(1 - a)Xip2 + aX2pi + ^(c/3i - bf32)Bi2 - . (71) 

Additional linear equations for the integral variables can be obtained from the previous ones and the 
following ones for the derivatives, 

H' = -{^B,+uB,2 + ^B2), (72) 

X[ = -(uB^+bBu), (73) 
X^ = -{CB12+UB2). (74) 

These follow taking the derivative of l|55|l . (|59|1 and H6Q(I with respect to L2. Therefore all the integral 
variables can be written in terms of the functions u, b and c by solving a linear system. In particular we 
get 

H = {Wna {1 - a) m)~Uu^ . (75) 
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2.4 Equations for the partition function and identification with Sine-Gordon 
model correlators 

Combining and l(75|l we have 



LOl log Za 



dy y u^{y) ■ 



(76) 



Comparing (|52() and H76|l with the results of we see that the partition function satisfies the same 
equations than the inverse of a correlator of exponential operators in the SG model at the free fermion 
point. In fact it is 

" ' ^^^^ 

where the correlator in the right hand side is evaluated in the Sine-Gordon model with action 



A 



(fx 



(78) 



and Va is the operator : exp[i\/47ra0] :. From this identification we learn that the long distance behaviour 
of u is 

u{t)~^- sm{an)Ki^2a (t) ■ (79) 

TT 

A correspondence between correlators of disorder operators for a scalar field with the inverse of 
correlators of exponential operators in the SG model was found in jllj . by exploiting the equations 
for the form factors. Our results, obtained by a different method, agree with that calculation, when 
identifying our partition function Za with the correlator of boson disorder operators with charges a and 
I -a. 



3 The entanglement and alpha entropies for free fields 

Summarizing the results, the final expression for the a-entropies for integer a = n is 

^ n — 1 

Cn = 2(1 -n) ^ ' 



k=l 



Wa = = dy y ul(y) , (81) 



with 



2 

Uaif) — sm{aTT)Ki-2a{t) as t ^ oo . (83) 

TT 

We have introduced the 1/2 factor in H8()(l because we want to present the results for a real scalar instead 
of a complex one. The function c„ is plotted in Figure 1 for some values of n. 

In the conformal limit i — s- the function u diverges and has a series solution of the form 

^ + - l)*log(i6) + Oit) . (84) 

tlog{tb) 

The full series is determined by the value of the integration constant b (which is a function of a). However, 
this is not known at present (except for the case a = 1/2 JHI)- From eq. (|5^^|l we get for the leading 
terms 

Wait) ^ -2a (1 - a) - + O {\og-\t)) , (85) 
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Figure 1: The c„ function for a real scalar field with, from top to bottom, n = 2, 3 and 50. At the origin 
they take the value ^4;^ and they decay exponentially fast for large t. 



and by analytic continuation of the sum in Ij80() 



cit) = \^^)^0{Xo,-\t)). (87) 

The first term corresponds to the conformal case for which there is a general result identifying c = C/3, 
where C is the Virasoro central charge Here it is C = 1 since we are considering a real scalar field. 

The leading term for large distances of the entropic c-function can be obtained from H83|) using the 
trick shown in P|. We get 

c{t) ^ -t Ki{2t) . (88) 

The analogous results for a Dirac field have been presented in a previous paper |5] . We copy here the 
corresponding formulae to facilitate the comparison with the bosonic case 

^ (n-l)/2 

Cn = -fT—n] ^ Wfc/„(i), (89) 

/oo 
dyyul, (90) 



__=_«)^ + „,.(i_„.) + __^, pi, 

2 

Ua —* — siii{Tra)K2a{t) as t — > oo . (92) 

TT 

The similarities are striking. The reason is that here there is also an identification with SG correlators 
given by 

Za ~ (VaV-a) ■ (93) 

The short and long distance behaviour of the Dirac c-function are 

c{t) ^ ^-^flog'{t) + 0{flog{t)) for t^O, (94) 

c(t) - ^tKi{2t) for t^oo. (95) 

We see that the bosonic c-function interpolates between the ones corresponding to a Dirac field (which 
has the same central charge at the conformal point) and a Majorana one (which has half the Dirac field 
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Figure 2: The entropic c-function c{t) = LdS{L) / dL for a real scalar field (dotted line) and the Dirac and 
Majorana fermion ones (top and bottom dashed lines respectively). The bosonic c-function interpolates 
between the Dirac one at the origin, where both tend to the conformal value c = 1/3, and the Majorana 
c-function for large i, where they decay exponentially fast. The cusp at the origin of the bosonic function 
is due to a l/log(i) term. These curves were obtained by numerical computation on a lattice of up to 
600 points. 

entropy). However, it is evident that a deeper understanding of this phenomenon requires more than the 
explicit calculation presented here. We have plotted the entropic c-functions for bosons and fermions in 
Figure 2. This curves are calculated numerically putting the model in a lattice and using the method 
given in Appendix A. Note also the very different behaviour at the origin. The bosonic c-function has a 
l/log(i) term that can be ascribed to the zero mode which is present at the conformal point. In fact, for 
< — > this behaviour gives for the entropy 

S^^ log(L/e) + i log (- log(m e)) ~ ^ log (- log(m L)) , (96) 

where e is an ultraviolet cutoff. In addition to the well-known logarithmic ultraviolet divergence, this last 
expression shows an explicit infrared divergence. 

4 Universal terms for the entanglement entropy in arbitrary 
dimensions 

Some universal terms for the entanglement entropy in more than two dimensions can be obtained from the 
two dimensional results. Specifically, let us consider a d-dimensional box, where d is the space dimension. 
Take one side of length L lying along the Xd axes direction to be much shorter than all the other sides 
having lengths i = l,...,d— 1, and lying parallel to the Xi axes directions. The entanglement entropy 
of the box is divergent with a term proportional to the (d — l)*'^-power of the cutoff and the total surface 
area E] ■ Several other divergent terms proportional to lower dimensional geometric quantities also 
appear in the entropy, including a term with universal coefficient, proportional to the logarithm of the 
cutoff, which is due to the sharp angles at the box vertices However, one expects to find also 

subleading contributions which are finite and universal. In particular in the limit Ri — > oo one expects a 
leading universal term having the form 

S ~ non-univ — ^ j^^-i ' ^^'^^ 

where the area A = IliRi, and fc is a universal and dimcnsionlcss function of L and the renormalized 
parameters of the theory alone. For free bosons and fermions we can calculate exactly the function k. 
In the limit of large Ri we can compactify the Xi directions, with i — 1, ...,d — 1, by imposing periodic 
boundary conditions Xi — xi -f Ri. As the term we are looking for is an extensive quantity in these 
coordinates, the compactification cannot change the result. For a real scalar field of mass m we Fourier 
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decompose it into the corresponding field modes in the compact directions. The problem then reduces 
to a one dimensional one with an infinity tower of massive scalar fields 4>qj^,...,q^_i, where qi,...,qd-i are 
integers. These have masses given by 

=-^ + E(|'Z.)'. (98) 

Thus, summing over all the decoupled one dimensional fields, one has 

dS{L) 



dL 



^c(LM,„...,,,_J , (99) 



with c{t) the one dimensional bosonic c-function. As we consider the large Ri limit we convert the sums 
into integrals and obtain 

m - ^^^^ t'-' dy, y^' £ dy, yt' c (^^J yf + yj^ , (100) 

where vo^iS"*^^) is the volume of the [d— 2) dimensional unit sphere (vo^iS*^) = 2). The same expression 
holds for free fermions. In this case the one dimensional c-function corresponds to the one of a Majorana 
field, and there is an additional factor of the dimension of the spinor space. 
In particular, in the massless limit and in 2 + 1 dimensions we have 

1 r°° 

ks ^ - dtcs{t) - 0.039, (101) 







1 f°° 

ko = - dtcoit) - 0.072. 
Jo 



(102) 



for a real scalar (S) and a Dirac fermion (D) respectively. The values on the right hand side were found 
by numerical integration of the curves of Figure 2. We have checked these last results numerically on a 
two dimensional lattice for rectangular sets of size up to 10 x 100 points, finding a perfect accord within a 
3% error. Further details of this and other results on universal terms for the entropy in 2 + 1 dimensions 
will be published elsewhere |14j . 



5 Conclusions 

In this work we have completed the discussion on the alpha and entanglement entropies for free massive 
fields in two dimensions, presenting the results corresponding to a scalar field. The Dirac fermion case 
was analyzed in a previous paper by a different approach. Here we arrive to the exact expression for 
the a-entropies for integer a and expansions for short and long distances of the geometric entropy for a 
massive scalar field. We do it by direct calculation of the relevant information on the Green function on a 
plane, where boundary conditions are imposed on a finite interval. In both cases, bosonic and fermionic, 
we find a relation to correlators of exponential operators in the Sine-Gordon model, associated to the 
same type of non linear differential equation of the Painleve V type. 

As we previously noted for the fermionic case '5', also for the scalar field the expression for the Cq 
functions show that they are decreasing in an explicit manner (see (|81|) and (|90|l ). This could be a hint 
to a possible extension of the entanglement entropy c-theorem. 

The method used in this work seems to be very well suited to generalizations, and thus it might be 
useful to calculate different contributions to the geometric entropy in more dimensions. It essentially gives 
information from the translation and rotation symmetries of the wave equation, even in the presence of 
symmetry breaking boundary conditions, once one has understood the structure of the singularities. 

Concerning the geometric entropy in arbitrary dimensions, we have also discussed universal terms 
which are related to the two dimensional results. We have found a universal term of the form 

Su - -fc-A- (103) 
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for the entanglement entropy of a thin box in d spatial dimensions, with a small side of length L and 
big phase area A. In some sense, this term could be interpreted as a universal (cutoff independent) form 
of the area law for the entropy, with L playing the role of the cutoff. It is known 7 that in quantum 
field theory a universal measure of information must be expressed in terms of the mutual information 
I{X, Y) ~ S{X) + S{Y) — S{X U Y) between two different non intersecting sets X and Y, rather than 
the entropy for a single set. In this sense it is remarkable that the term (|103ll determines the mutual 
information between to big boxes as their parallel phases come into contact. 
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7 Appendix A: numerical entropy in a lattice 

We use the method presented in ^3] to give an expression for pA in terms of correlators for free boson 
and fermion discrete systems. The case of a fermion field is summarized in 0. For the scalar field it is a 
generalization of the method given by Sredniki . We start with a free Hamiltonian for bosonic degrees 
of freedom with the form 

where (/j^and tt^ obey the canonical commutation relations [0i,7rj] = iSij^ and M is a Hermitian positive 
definite matrix. The vacuum (ground state) correlators are given by 

X, = (<^,;0,) = i(M-^)y, (105) 
R^ - (^,7r,) = i(M^),,. (106) 

Let X^ and -Pj^be the correlator matrices restricted to the region A, that is i,j G A. Due to the Wick 
theorem any other correlator inside A can be written in terms of these two. The expectation values 
calculated with pA must also satisfy Wick's theorem and give place to the two point correlators X^ and 
P,^. This fixes pA to be of the form 

PA = Ke^^""'"' , (107) 

where K is a normalization constant and the a; are independent boson annihilation operators, [ui, a|] = 
Sij, which can be expressed as linear combination of the and tTj, i, j G A. The ei are related to the 
eigenvalues i/i of \/X^.P^ by 

tanh(|)^i-. (108) 
Hence, the entropy can be evaluated as Sa — ^ Si, where 

Si = -log(l -e-O +e;Y3^ = + + ^) - - - ^) ■ (109) 

Thus, to compute the entropy numerically we need to diagonalize the matrix X^.P'^. Its eigenvalues 
lyf, due to the uncertainty relations, are always greater than 1/4. Note that for the total system XP = j 
and vi = ^, with zero entropy. 

The lattice Hamiltonian for a real massive scalar is given by 



1 



n=0 
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We have taken TV lattice sites and set the lattice spacing to one. The correlators 1105() and (|106|l for the 
infinite lattice limit N ^ oo are 



dx- 



2^m2 + 2(1 - cos(27rx)) ' 
1 1 

2.^x(rn-n) ^^2 + 2(1 - C0s(2^x)) . 



(Ill) 

(112) 
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